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We determine an upper bound for a possible right-handed 6— >c quark current admixture in semilep- 
tonic B — >■ Xci'i' decays from a simulateous fit to moments of the lepton-enery and hadronic-mass 
distribution measured as a function of the lower limit on the lepton energy, using data measured 
by the BaBAR detector. The right-handed admixture is parametrized by a new parameter cji as 
coefficient of computed moments with right-handed quark current. For the standard model part 
we use the prediction of the heavy-quark expansion (HQE) up to order 1/mf and perturbative cor- 
rections and for the right-handed contribution only up to order 1/ml and perturbative corrections. 
We find = 0.05 ^p gp in agreement with the standard-model prediction of zero. Additionally, 
we give a contraint on a possible right-handed admixture from exclusive decays, which is with a 
value of c'n = 0.01±0.03 more restrictive than our value from the inclusive fit. The difference in 
|Vcbl between the inclusive and exclusive extraction is only slightly reduced when allowing for a 
right-handed admixture in the range of c'n — 0.01±0.03. 



I. INTRODUCTION 

Parity violation is implemented in the Standard Model 
(SM) by assigning different weak quantum numbers to 
left- and right-handed quarks and leptons and it is fair 
to say that there is yet no deeper understanding of the 
symmetry breaking mechanism in weak interactions with 
respect to parity transformations. On the experimental 
side, parity violation is well established in the leptonic 
sector, e. g. through the measurements of the Michel pa- 
rameters in the decay of Muons. 

However, in hadronic transitions parity violation is 
much harder to test due to the uncertainties present in 
the calculation of the hadronic matrix elements of the 
quark currents. In turn, this leaves sizable room for a 
possible non-left-handed admixture. While this is mainly 
true in the case of light quark systems, the calculational 
methods have significantly developed for heavy quarks 
using the fact that the heavy quark masses are large com- 
pared to the binding energy of heavy hadrons. 

In particular for inclusive semileptonic &— )-c transitions 
the theoretical methods are in a mature state and are 
applied in the framework of the SM to extract e. g. the 
CKM Matrix element |Vcfc| with an unprecedented rela- 
tive precision of less than 2 % [Ij . 

Clearly the precision of the data as well as of the the- 
oretical methods may serve also to perform a test for 
non-standard couplings. In two recent papers [51 [3] the 
necessary calculations have been performed to check for 
a non-left-handed coupling in inclusive semileptonic &— >-c 
transitions. In the present paper we use these results to- 
gether with the BaBAR data to obtain information on a 
possible right-handed admixture to the weak &— >-c cur- 
rent. 

Recently the tension between the inclusive and exclu- 
sive determinations of \ Vub\ and — to a lesser extend — also 
of \Vcb\ motivated speculations to explain this by right- 



handed admixtures in the weak hadronic currents. In 
[4] it is shown that a right-handed admixture can soften 
the tension and that a right-handed admixture can be 
obtained within the MSSM. 

In the next section we recapitulate the theoretical in- 
put and fix our notation. In section III we perform the 



analysis based on BaBAR data. In section |IV| we discuss 
our result and compare them in section |V] to limits from 
exclusive decays. Finally, in section IVl] we conclude. 



II. THEORY BACKGROUND AND NOTATION 

It is well known that any new physics effect beyond the 
SM can be parametrized in terms of higher- dimensional 
operators, which are singlets under the SM symmetry 
SU(3)QCDxSU(2)wcak X U(1)y. Assuming that the Higgs 
sector is minimal (i.e. if one considers only a single Higgs 
doublet) there is only one operator at dimension five 
which is related to a Majorana mass of the neutrino and 
hence only affects the leptonic sector. At dimension six 
one finds a long list of possible operators i5| among which 
we find also operators modifying the helicity structure 
that appears in the semileptonic 5— >c decays. 

After spontaneous symmetry breaking SU(3)qcd x 
SU(2)woak X U(1)y SU(3)qcd X U(l)cm and after run- 
ning down to the scale of the 5-quark mass one finds for 
the effective interaction [31 [3] 
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where Jj^ — ej'^P- Ve is the usual leptonic current and 
J\^,^l is the generalized hadronic &— >c current which is 
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given by 

+ g-L cit)lp^b + g^ cii)lp+b 

+ di^ id^{cia^t,P-b) + id''{cia^i,P+b) , 



(2) 



where P± denotes the projector on positive/negative chi- 
rahty and is the QCD covariant derivative. Note that 
the leading term contributing to the rate wiU be the in- 
terference term with the SM (cx clCr), which means that 
the leptonic current remains as in the SM since we con- 
sider only final states with electrons and muons and thus 
can neglect the lepton mass, 

Furthermore, cl contains the SM contribution and 
hence cl = l + 0{v'^ / A^), where v is the vacuum expecta- 
tion value from spontaneous symmetry breaking and A is 
the new-physics' scale. All other contributions can only 
appear through a new-physics effect. In particular, the ef- 
fective field theory approach reveals that cr = 0{v'^/A'^), 
while all helicity changing contributions are expected to 
be further suppressed by a small Yukawa coupling [2 E] . 

In the following analysis we restrict ourselves to an in- 
vestigation of the parameter cr. As has been shown in [3] 
the lepton-energy moments and hadronic-mass moments 
are not very sensitive to the parameters ^l, 5r, (^l and 
dR. Because the moments depend on the squared ma- 
trix element the parameters appear in pairs, of which the 
leading contributions are cf^ and clCr. For the combined 
fit the parameter cl can also be dropped as an overall 
factor being absorbed in \Vcb\- Thus the parameter used 
in the fit is Cr = cr/cl. 



III. ANALYSIS 

A. Fit Setup 

The combined fit for the extraction of the new param- 
eter Cr is performed along the lines as described in [7] 
using the HQEFitter package |H]- It is based on the 
minimization, 
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with the included measured moments Afoxp, the corre- 
sponding theoretical prediction of these moments Mthco 
and the total covariance matrix Ctot defined as the sum 
of the experimental ( Ccxp) and the theoretical ( Cthco) 
covariance matrix, respectively. 

In the analysis of 7] the theoretical prediction for 
the moments A/hqe are calculated perturbatively in 
a Heavy-Quark Expansion (HQE) in the kinetic-mass 
scheme up to 0{l/ml) with perturbative contributions 
[SUllj resulting in a dependence on six parameters: the 
running masses of the b- and c-quarks, mf,(/i) and TOc(^), 
the parameters fi^ and at 0{l/ml) in the HQE, and, 
at 0{l/ml), the parameters and pj^g. 



New in this analysis is the inclusion of possible right- 
handed quark currents in the calculation of the theoret- 
ical prediction of the moments. The right-handed con- 
tributions are calculated and used here up to 0{l/ml) 
in the HQE and 0{as) in the perturbative correction. 
The aim of this fit is to give an upper bound for the 
relative contribution of a right-handed current compared 
with the standard-model left-handed current, which is 
parametrized by a prefactor Cr for the new contributions 
to test. Thus the theoretical prediction of the moments 
depends on seven parameters to fit: 

Mtheo = Mtheo (Cr , TOb , mc , /X^ , , pI, , pis ) . 



B. Determination of |V^b| 

In the presence of a right-handed mixture the definition 
of the parameter \Vcb\ becomes ambiguous. Out of the 
three parameters \Vcb\, ci,, Cr only two are independent, 
since cl can be absorbed into \ Vcb\- To this end we choose 
to define 



\Vcb\bLlp.CL 



\Vcb\ {bLlfiCL + d^bR-f^cii) 



For the determination of |Vcb| the fit uses a linearized 
form of the semileptonic rate T sl expanded around a- 
priori estimates of the HQE parameters [S]: 



Bciv 1.55 



0.0417 V 0.1032 + ^.30 (a^b) - 0.22)] (4) 

X [l - 0.66 (mf,- 4.60) 0.39 (mc- 1.15) 
+ 0.013 (pI ~ 0.40) + 0.09 (pI, ~ 0.20) 
+ 0.05 {p% - 0.35) - 0.01 {pis + 0.15) 
+ 0.341 4] . 

Note the last term (0.341 Cr), taking into account the 
possible contributions from a right-handed quark current. 
The a-priori estimate of Cr is zero, i. e. the standard- 
model value. Due to the sizable factor and positive sign, 
a positive value of Cr increases \Vcb\ compared to the 
standard-model fit without c^. 

The total branching fraction B{B — Xci^V) in the fit 
is extrapolated from measured partial branching fractions 
Bp'^^.^{B -> Xct~V), with pI > p'l^in- This is done by 
comparison with the HQE prediction of the relative decay 
fraction (r.h.s): 
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Thus the total branching fraction can be introduced as 
a free parameter in the fit. By adding the average B- 
meson lifetime tb (average between neutral and charged 
_B-mesons, see also next paragraph) to the measured and 
predicted values, Vcb\ can as well be introduced as a free 
parameter using (4]). 
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Figure 1: The measured lepton-moments (•/o) compared with the result of the simultaneous fit (solid red line) as function of 
the minimal lepton momentum pJmin- The measurements included in the fit are marked by solid data points (•). The dashed 
lines indicate the theoretical fit uncertainty obtained by the variation of the fit parameters in order to convert their theoretical 
uncertainty into an error of the moments. 




Figure 2: The measured hadronic-mass moments (•/o) compared with the result of the simultaneous fit (solid red line) as 
function of the minimal lepton momentum p^^min- The measurements included in the fit are marked by solid data points (•). 
The dashed lines indicate the theoretical fit uncertainty obtained by the variation of the fit parameters in order to convert their 
theoretical uncertainty into an error of the moments. 



C. Experimental Input 

The combined fit is performed vifith a selection of the 
following 25 moment measurements by BaBAR which are 
characterized by correlations below 95 % to ensure the 
invertibility of the covariance matrix: 

• Lepton energy moments measured by BaBAR 
[T^ . We use the partial branching frac- 
tion Sp« . at the minimal lepton momentum 

> 0.6, 1.0, 1.5GeV/c, the moments (Ef,) for 

> 0.6,0.8,1.0, 1.2, 1.5 GeV/c, the central mo- 
ments {{Ei-{Ei,)f) for p| > 0.6, 1.0,1.5GeV/c 
and {{Et-{Ei)f) for p| > 0.8, 1.2 GeV/c. 

• Hadronic mass moments measured by BABAR 
[7]. We use the moment {m\) for 

> 0.9, 1.1, 1.3, 1.5 GeV/c and the central mo- 
ments ((to^ ~ ('^x))^) (("^x ~ ('^x))^) 
both for p| > 0.8,1.0, 1.2, 1.4 GeV/c. 

Furthermore we use the average B meson lifetime tb — 
foTo + {l-fo)r± = (1.585±0.007) ps with the lifetimes of 
neutral and charged B mesons tq and t± and the relative 
production rate, fo = 0.491 ± 0.007, as quoted in [T5] . 



D. Theoretical Uncertainties 

Theoretical uncertainties for the prediction of the mo- 
ments Mthoo are estimated by variation of the parame- 
ters. The standard model parameters, that are all ex- 
cept c^, are treated as in [7]. The uncertainty in the 
non-perturbative part are estimated by varying the cor- 
responding parameters /ij and by 20 % and and 
Pls tiy 30 % around their expected value. For the uncer- 
tainties of the perturbative corrections — 0.22 is var- 
ied up and down by 0.1 for the hadronic mass moments 
and 0.04 for the lepton energy moments and the uncer- 
tainties of the perturbative correction of the quark masses 
mil and rric are estimated by varying them 20 MeV/c'^ up 
and down. An additional error of 1.4% is added to \Vcb\ 
from the fit for the uncertainty in the expansion of the 
semileptonic rate Fsl, which is not included in the fit, 
but quoted separately as theoretical uncertainty on | Vcb \ ■ 

Additionally the influence of the right-handed contri- 
butions on the theoretical uncertainties in the predictions 
of the moments has to be included. Varying in a sim- 
ilar fashion as the other parameters, around the a priory 
estimate of zero showed only very little influence on the 
fit results. Due to the fact that the right-handed contri- 
butions are included up to 1/to^ in the non-perturbative 
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and 0{as) in the perturbative corrections for all mo- 
ments, the uncertainties in the prediction of the moments 
are not sizable and thus the variation of has to be 
rather small. For the final results the variation of has 
not been included, because of no influence on the signif- 
icant digits. 



IV. RESULTS 

Table U shows the fit results and table HIl the corre- 
sponding standard-model fit results, which were obtained 
by performing the fit with fixed to zero. Figs, [l] 
and [2] show a comparison of the fit results with the mea- 
sured moments for the lepton-moments and the hadronic- 
mass moments, respectively. The uncertainties Aoxp and 
Athcor are the expected experimental and theory errors 
determined by Toy-MC studies (see [T3]) while Atot is 
the total uncertainty provided by the fit. 

The estimate for = 0.05 ^[j'^p is consistent with the 
standard-model prediction of zero, but the uncertainty 
reveals an unexpected low sensitivity of the semileptonic 
fit to possible right-handed contributions. We state the 
upper relative admixture limit of |cp| = 0.9 at 95% con- 
fidence level. 

The extracted value of \Vcb\ = (43 +g) • 10"^ is con- 
sistent with the value from the standard-model fit, but 
its uncertainty is quite difi^erent. In our fit, this is due 
to the influence of a sizable Cp uncertainty on the de- 
termination of I Vcb I in ([4 ) , which becomes evident in the 
contour plot of the (c^Vcbl) plane, showing a shallow 
and steep covariance ellipse. 

To compare the quality of the fits the P-value 
(prob(x^, ''^dof)) suits best, because the fits differ by 
their number of degrees of freedom and thus the 
value alone is not sufficient. For the fit with we 
find = 7.299 with 17 degrees of freedom and thus 
prob(7.299, 17) = 0.979 and for the standard-model fit 
prob(7.312, 18) — 0.987, which shows neither improve- 
ment nor worsening. 

The uncertainty of the result for is dominated 
by the theory error Athoo=0.42 in table |l] compared to 
Acxp=0.13. As a consequence including additional exper- 
imental data, e. g. from Belle, will not improve the limit 
for a possible right-handed contribution at this point. We 
investigated the theory error by changing the variation of 
the parameters as described in section [TlI D[ It turned out 
that the theory error cannot be pinned down to the uncer- 
tainty of a specific parameter. Furthermore, decreasing 
all theory errors had only little effect on the theory error 
of c'^ . We come to the conclusion that the shape of the 
considered spectra and hence their moments are too sim- 
ilar for a left and right-handed 6 — > c current, ending up 
in a low sensitivity of and a weak constraint therein. 
Thus, improving the theoretical description, either for 
the standard-model part or the right-handed contribu- 
tions, e.g. including the l/m^ and BLM corrections of 
the right-handed current, will not reduce the uncertainty 
of Cr significantly in the performed moment analysis. 



V. RIGHT-HANDED ADMIXTURE FROM 
EXCLUSIVE DECAYS 

It is interesting to note that the value of \Vcb\ ex- 
tracted in this way \Vcb\ = (43 ^g) • 10~^ is in agree- 
ment with I Vcb I from exclusive decays, from which 
\Vcb\ = (38.6±1.3) • 10-3 (from B DH-v) 15J is 
obtained, while the value from the standard-model fit 
\Vcb\ = (41.9±0.8) • 10-3 is not. This is due to the low 
sensitivity of and thus the large uncertainty of the 
extracted value. In addition, also the exclusive decays 
allow us to constrain a possible right-handed admixture 

Hi HI]. 

The most straightforward way of obtaining this infor- 
mation is to study the exclusive differential rates at the 
point of maximal momentum transfer to the leptons, cor- 
responding to equal four-velocities of the initial and fi- 
nal hadron. We consider the decays B — )• D£~i' and 
B — )• D*l^v. The corresponding rates in the standard 
model are usually parametrized in terms of two form fac- 
tors; the relevant expressions close to the point of maxi- 
mal momentum transfer read 



dr 



B^D 



Aw 

dr^^ 

dw 



= rol6r3(r + lf{w^ - lf'^mb\G{w)) 
= rol92r3(r, - lf{w^ - if'WV^bm^)) 



2 

(6) 



where w=v-v' is the scalar product of the hadronic veloci- 
ties, r=raD /mB, rt,=mo* Ims, and ro=G|.m^/(1927r3). 
The information which is extracted by the experiments 



Vcb\ 



lim 



dr 



B^D 



1 



lim 



w^i dw rol6r3(r+ l)2(u)2 - 1)3/2' 

dr^^^* 1 



dw rol92r3(n ~ l)2(u;2 - 1)1/2 



(7) 



which in the standard model is the product of the form 
factors at It; = 1 and \ Vcb\- Combining this with a theoret- 
ical prediction of the form factors ai w = 1 one extracts 
\V.b\- 

At the non-recoil point w = 1 the B ^ D transition 
is completely dominated by the vector current, while the 
B — ^ D* decay is proportional to the axial vector current. 
Thus, including a right-handed admixture, the informa- 
tion extracted from ([7| is |cl + cr||Vc6|^(1) for the case of 
the B ^ D transition and |cr — cl| | Vcb] J^(l) for B ^ D* . 
The current experimental data yield |18| : 

\ci^ + CT^\\Vcb\Q{l) = (42.4 ±1.56) X 10-3 (8) 
lcR-CLllKhl-^(l) = (35.41 ±0.52) X 10-3 (9) 
Using the lattice data (which are also used to extract 

\Vcb\) [MS] 



g{l) = 1.074 ±0.024 
=0.921 ±0.025 



(10) 
(11) 
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Table I: Results of the full fit for c'^ and the canonical set of parameters \Vcb\, rnt, rric, B, fi^, fj,c, Pd B,nd p\g, separated by 
experimental and theoretical uncertainties. For IVcbl we take into account an additional error of 1.4% for the uncertainty in 
the expansion of the semileptonic rate Vsl- Correlation coefficients for the parameters are listed below. The uncertainties Aexp 
and Athcor are the expected experimental and theory errors determined by Toy-MC studies (see [TJ) while Atot is the total 
uncertainty provided by the fit. 









mi, 


rric 


B 


2 


2 


Pl 


Pls 






xlO~^ 


[GeV/c^] 


[GeV/c^] 


[%] 


[GeV^] 


[GeV^] 


[GcV^] 


[GeV''] 


Results 


0.0517 


42.61 


4.588 


1.133 


10.674 


0.472 


0.303 


0.192 


-0.122 


Aexp 


0.1335 


2.00 


0.000 


0.000 


0.241 


0.032 


0.052 


0.015 


0.096 


Atheo 


0.4209 


5.55 


0.110 


0.149 


0.068 


0.094 


0.035 


0.056 


0.012 


Arsi 




0.60 
















Atot 


+0.3356 


+4.76 


+0.081 


+0.112 


+0.240 


+0.121 


+0.061 


+0.095 


+0.095 




-0.4962 


-6.34 


-0.158 


-0.226 


-0.274 


-0.086 


-0.064 


-0.043 


-0.099 


4 


1.00 


0.99 


0.78 


0.74 


0.67 


-0.79 


0.33 


-0.84 


0.21 






1.00 


0.75 


0.72 


0.72 


-0.77 


0.29 


-0.82 


0.22 


nib 






1.00 


0.99 


0.64 


-0.72 


0.28 


-0.68 


0.17 


nic 








1.00 


0.65 


-0.70 


0.16 


-0.64 


0.20 


B 










1.00 


-0.47 


0.16 


-0.52 


0.14 


2 












1.00 


-0.21 


0.84 


-0.14 


Pa 














1.00 


-0.31 


0.03 


P% 
















1.00 


-0.29 


pis 


















1.00 





Table 11 


Results of the standard model fit with 


c^Yi fixed to 


zero. 








\V.b\ 


nib 




B 




Mg 


Pl 


pls 




xlO"^ 1 


GeV/c^] 


[GeV/c^] 


[%] 


[GeV^] 


[GeV^] 


[GeV^^] 


[GeV''] 


Results 


41.93 


4.578 


1.120 


10.654 


0.482 


0.300 


0.198 


-0.125 


Aexp 


0.44 


0.058 


0.085 


0.175 


0.023 


0.040 


0.015 


0.081 


Atheo 


0.38 


0.045 


0.067 


0.061 


0.055 


0.043 


0.027 


0.049 


Ar,, 


0.59 
















Atot 


+0.83 


+0.074 


+0.106 


+0.185 


+0.058 


+0.059 


+0.031 


+0.094 




-0.83 


-0.070 


-0.107 


-0.185 


-0.060 


-0.059 


-0.031 


-0.094 


\v.A 


1.00 


-0.35 


-0.21 


0.67 


0.29 


-0.39 


0.35 


0.06 


rub 




1.00 


0.98 


0.25 


-0.26 


0.06 


-0.08 


0.01 


rric 






1.00 


0.29 


-0.29 


-0.11 


-0.07 


0.07 


B 








1.00 


0.14 


-0.08 


0.10 


0.00 


pl 










1.00 


0.08 


0.52 


0.04 


2 

Pa 












1.00 


-0.07 


-0.04 


Pl 














1.00 


-0.20 


pls 
















1.00 




black contour in the (cj^, |Vc6|) plot shows the result from exclusive decays (cp — 0.01 ± 0.03) as computed in section [v| 
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we can extract the ratio Cp = cr/cl to be 
4 = 0.01 ± 0.03 



(12) 



with the assumption of no sizable correlations between 
the experimental measurements of the right-hand sides 
of Eqns. Q and ^ as well as between the form 
factor values given in ([To]) and ([IT]). The value for 
\Vcb\ extracted from Eqns. (|8| and ^ is found to be 
IKhloxci = (39. 0]*]^ ^) • 10^"^ which has to be compared 
to IKbloxci = (38.8 ± 1.0) • IQ-^ when setting ^ = in 
Eqns. (|| and (|9]). 

The result of is compatible with zero and, in fact, 
more restrictive than the determination from inclusive 
decays. Obviously the exclusive decay gives access to 
data separated by the handedness of the 6 — > c current in 
contrast to the inclusive decay, leading to a better limit 
on possible right-handed contributions. 

In turn, we can use the result for 4 to de- 
termine iVcblinci and compare with iVcfilcxci- This 
can be done by imposing a Gaussian constraint 
of Cr=0.01±0.03 in the fit with a possible right- 
handed current. The result |K6|inci=(42.0±0.9)-10-3 
compared to |Vcb|cxci=(39.0^Jg)-10^'^ exhibits a ten- 
sion by 3.0- lO^'^ of the central values. Determin- 
ing \Vcb\ by ([S]) and ^ with Cr set to zero gives 
|V'ch|exci(cR=0) = (38.8±1.0)-10~^ and allows us to exam- 
ine the differences in the tensions between inclusive and 
exclusive decays with and without a right-handed cur- 
rent, by comparing this value to the standard model 
fit value |ycb|mci(cR=0)=(41.9±0.8)-10-3 (see Table |ll]) 
yielding a tension of about 3.1-10"^ of the central values. 
As a consequence, the difference in the central values 
between \Vcb\ exclusive and inclusive is slightly reduced, 
and more importantly, the uncertainty on |Vc&| exclusive 
is considerably larger when allowing for a right-handed 
admixture resulting in a smaller significance of the ob- 
served effect. In our analysis, which is using only the 
inclusive BaBAR data, the difference between exclusive 
and inclusive is reduced from a 2Aa to a 2.1a effect. 



VI. SUMMARY 

We have performed a full-fledged fit to moments of the 
lepton-energy and hadronic-mass distribution of semilep- 
tonic B — >■ Xc£~i' decays, including a possible right- 
handed admixture to the b ^ c current. We have con- 
sidered the non-standard contributions up to I /ml in 
the non-perturbative and 0{as) in the perturbative cor- 
rections. The corresponding fit in the framework of the 
standard-model yields the most precise determination of 
|Vc&|, due to the elaborated theoretical description and 
the precise measurements of the B factories [TS]. Our fit, 
including a right-handed admixture, is in agreement with 
the standard model assumption of zero for a right-handed 
contribution. Unfortunately, the result Cr — 0.05 ^pg^ 
reveals a low sensitivity of the fit to a right-handed con- 
tribution, compelling us to state the upper relative ad- 
mixture limit of |cr| = 0.9 at 95% confidence level. The 
moments of the spectra used in the fit are too similar 
for right- and left-handed contributions, resulting in the 
quoted low sensitivity and weak bound of Cr. 

Exclusive decays are competitive in the determination 
of I Vcb I , given the precise values for the form factors at the 
non-recoil point obtained from lattice QCD calculations. 
The same precise values from lattice QCD calculations 
can be used to obtain a constraint on Cr, which is con- 
siderably stronger than the one obtained from inclusive 
decays: Cr=0.01±0.03. Using this result to determine 
I V'cb I inci we can compare the tension between | Vcblinci and 
|V'cf,|cxci without a right-handed current and with a right- 
handed current contribution as from exclusive decays. A 
right-handed current reduces the tension by about 3 % 
and its significance from a 2.4 cr to a 2.1 cr effect. 
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